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The suppression of quarkonium production in nucleus-nucleus collisions was originally 
proposed as a signal of colour deconfinement. Strong "anomalous" J/ip suppression 
in Pb-Pb collisions has been reported at this Conference by the NA50 Collaboration. 
Is this suppression really anomalous? Can we conclude that the quark-gluon plasma is 
already discovered? What has to be done next? I address these questions basing on the 
(N| ■ current theoretical understanding of quarkonium production and new precise experimental 
information. 

\o : 

1. INTRODUCTION 

' 

Q_i! The idea to use heavy quarkonia as a probe of excited QCD matter produced in rel- 
ativistic heavy ion collisions was proposed a decade ago This suggestion was based 
on the concept of colour screening of static potential acting between the heavy quark 
and antiquark, which occurs in hot and/or dense quark-gluon matter. A year later, the 
J/tp suppression in nucleus- nucleus collisions was observed experimentally @. Since that 
time, quarkonium suppression has always been a respectable prospective signal of de- 
confinement. However, the experimentally observed suppression has been consistent not 
only with the deconfinement scenario, but with a plenty of "conventional" explanations 
as well [H]. Moreover, the differences between various conventional approaches, invoking 
such different mechanisms as nuclear absorption, gluon shadowing, or comover absorption, 
but nevertheless all providing more or less reasonable description of the data, created a 
lively controversy. The problem thus was clearly awaiting a more detailed and system- 
atic analysis, many essential inputs for which were missing. In fact, the lack of precision 
data on J/ijj and ip' production in p — A collisions, poor theoretical understanding of J/ ip 
production and its interactions with light hadrons made this analysis virtually impossible. 

Fortunately, the situation has started to change recently: new high precision data on 
quarkonium production (see and references therein) have significantly advanced 

the theory J7j , and the operator product expansion techniques have allowed a systematic 
calculation of quarkonium absorption cross sections fj-jlOll and dissociation rates in con- 



fined matter |]I0||- ||14|| . The dominance of higher Fock states in quarkonium production, 
revealed by the new Fermilab data || and naturally emerging theoretically, has inspired 
a new approach to nuclear attenuation of quarkonium production ||15| . This approach, 
as I shall discuss in this talk, has enabled a quantitative understanding of quarkonium 
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suppression in both p — A and A — B collisions, giving a credit to the nuclear absorp- 



tion model |IB |. It has become clear that the existing S — U data in fact do not exhibit 
any anomalous behaviour and are not consistent with a deconfinement scenario, which 
requires additional strong suppression. This is why the results on J/ip suppression in 
Pb-Pb collisions were so anxiously awaited - they were our last chance to see something 
unusual before the advent of future experiments at RHIC and the LHC 

These results, presented by the NA50 Collaboration at this Conference |17| , are striking. 
The data clearly show a strong J / tp suppression, going way beyond the expected. But is 
this suppression really anomalous? Have we finally reached the border of the long-awaited 
terra incognita of deconfined quark-gluon matter? In this talk, I shall attempt to address 
these questions. 



2. QCD ATOMS IN EXTERNAL FIELDS 

2.1. Quarkonium Interactions and the Operator Product Expansion 

In the Operator Product Expansion (OPE) approach, the amplitude of heavy quarko- 
nium interaction with light hadrons is represented in the form 



F* h = if d*xe i<lx (h\T{J(x)J(0)}\h) = 5> n (Q,m Q )(0„>, 



(1) 



where the set {O n } should include all local gauge-invariant operators expressible in terms 
of gluon fields; the matrix elements (O n ) are taken between the initial and final light- 
hadron states. The coefficients c n are expected to be computable perturbatively and are 
process-independent . 



Figure 1. A sample diagram describing quarko- 
nium interaction with a light hadron in the OPE 
scheme; dashed lines are the gluon propagators, 
ovals represent the quarkonium wave function, and 
the blob stands for the gluon structure function of 
the hadron. 



The Wilson coefficients c n were computed for S and P |L4| , |18j states in the leading 
order in 1/N 2 (N is the number of colours). The expectation values (O n ) of the operators 
composed of gluon fields can be expressed as Mellin transforms [0 of the gluon structure 
function of the light hadron, evaluated at the scale Q 2 - - ~ 



(On) 



dx x n g(x, Q = t\ 



(2) 



Since the total $ — h cross section is proportional to the imaginary part of the amplitude 
F<5>h, the dispersion integral over the c.m.s. energy A leads to the set of sum rules, relating 
the cross section to the gluon structure function of the light hadron. This relation, illus- 
trated in Fig. 1, has a very important property: the magnitude and energy dependence 
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Figure 2. J /if) photoproduction cross section; the curve is the theoretical prediction p0| . 



of the quarkonium dissociation cross section at low energies is entirely determined by the 
behaviour of the gluon structure function at large x ~ 1/A, whereas the cross section at 
high energy is governed by the small x behaviour of the structure function. Since the 
gluon structure functions of light hadrons are suppressed at large x, the calculated cross 
section rises very slowly from the threshold. When the hadron momentum in the J /if) 
rest frame is Ph — 5 GeV, the cross section is more than an order of magnitude below its 
asymptotic value. 

Recently, the calculation sketched above has been refined [p(J by taking into account 
target mass corrections, the real part of the scattering amplitude restored by dispersion 
relations, and the use of modern gluon structure functions inferred from the analyses of 
HERA data. This allows to evaluate the cross section in the entire energy range accessible 
to present experiments; the results confirm the threshold behaviour of the absorption 
cross section established previously. Vector meson dominance relates the cross sections 
of J /if) dissociation and photo-production; Fig. 2 shows the results compared to the 
available data. One can see that a strong threshold suppression of the J /if) absorption 
cross section is actually required by the data. 
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2.2. Quarkonium Interactions with Pions 

It can be shown that spontaneously broken chiral and scale symmetries of QCD imply 
decoupling of low-energy pions from heavy quarkonium fl4j. The proof is based on the 



application of low-energy QCD theorems |2T] (see ]22] for a recent review and introduc- 



tion) to the amplitude of quarkonium interactions with light hadrons [ 23fl . Qualitatively, 
the origin of the decoupling can be explained in the following way At low energies, the 
amplitude of quarkonium interaction is proportional to the gluon field operator dominat- 
ing the trace of the energy- momentum tensor of QCD. The appearance of this operator 
in the trace of the energy-momentum tensor is a reflection of the broken scale invari- 
ance of QCD, so the coupling is determined by the scale dimension of the hadron field. 
Chiral symmetry, however, implies zero scale dimension for the Goldstone boson fields - 
otherwise scale transformations would break chiral invariance. 

2.3. Quarkonium Production in Hadron Collisions 



The perturbative approach to quarkonium production ||24||,p5j| is based on the assump- 
tion that the production process is localized at distances ~ ■, much shorter than the 
size of quarkonium r ~ [a s (r~ l )mQ\~ l . This approach is justified if all gluons involved in 
the production carry a high momentum q ~ niq. However, the hadroproduction of vector 
states, for example, requires at least three gluons, of which only two must be hard to 
create the QQ pair. At small P? (the domain that dominates the integrated quarkonium 
production cross sections), the third gluon can be very soft, and is emitted (or absorbed) 
at distances of the order of quarkonium size. This is clearly inconsistent with the fac- 
torized form of the amplitude, and may "explain" the failure of perturbative approach 
in describing the integrated cross sections of quarkonium production at fixed target en- 
ergies. At collider energies, the perturbative approach fails even at high P?, since the 
non-perturbative contribution to the gluon fragmentation becomes important These 
arguments help to understand the phenomenological success of the colour evaporation 
model in explaining the data (see for a recent study). 

A consistent solution of this problem emerges if one assigns the soft gluon to the quarko- 
nium wave function introducing the notion of \QQg...) higher Fock states [/]]. In fact, such 
states appear naturally in the OPE scheme described above. Consider, for example, the 
amplitude of quarkonium interaction with an external gluon field (see Figs. 1 and 3): it 
includes the transformation of the colour- singlet quarkonium into a colour-octet 
state. The overall colour neutrality is of course preserved and ensured by the coloured 
gluon cloud surrounding the [QQ]s state. The simplest example of such system is provided 
by the \[QQ]sg) state. Since the vacuum of QCD has a complicated structure p7| with 
(g 2 G 2 ) 7^ 0, it induces a significant admixture of the IfQQJsfi 1 ) component in the wave 
function of quarkonium |28| - see Fig. 3a. For a physical J/ip state, this leads to the 



following generic decomposition: 

\J/ijj) = ai \cc) + a 2 \{ce] 8 g) + ... (3) 

Similar decompositions hold for other quarkonium states; for x states, for instance, the 
importance of higher Fock component is implied by the divergence of the perturbative 
annihilation amplitude in the soft gluon limit f29|. The magnitude of the \[QQ]%g) state 
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a) b) 

Figure 3. Interactions of quarkonium with external gluon fields; solid (dashed) lines are 
the heavy quark (gluon) propagators, and ovals represent the quarkonium wave function. 



admixture is reflected by the magnitude of relativistic corrections in the NRQCD approach 
[ p9| and by the size of power corrections in the QCD sum rule approach [p^j . These 
corrections are generally not very large, making applicable the familiar concept of heavy 
quarkonium as of a non-relativistic system essentially composed of just QQ state. However 
in certain processes - like production and annihilation of quarkonium - these components 
can play extremely important roleQ In fact, the leading order production of heavy vector 
quarkonium proceeds via the gluon fusion producing the QQ pair in a colour-octet state 
that later neutralizes its colour emitting (or absorbing) an extra gluon. If this extra gluon 
is soft (as is the case in the small Pt domain), the production process can be visualized 
as proceeding via the higher Fock state \[QQ]%g) (see Fig. 3b). 

Since the colour Coulomb interaction between the heavy quarks in the colour-octet state 
is repulsive and weak (~ 1/(N 2 — 1) with respect to the attraction in the colour-singlet 
state, where N is the number of colours), the ([(QQJsg) state is separated from the basic 
\QQ) state by the mass gap of ~ e , where e is quarkonium binding energy. This (virtual) 
state therefore has a proper lifetime of r ~ 1/eo- m the frame where quarkonium moves 
with momentum P, the superposition (|J) will be coherent over a distance z c ~ tP/2Mq. 
At high energies, this distance is sufficient for a produced \ [QQ]%g) state to traverse the 
entire nuclear volume. 

What will be the effect of the nuclear medium on the propagation of such a state? To 
answer this question, let us first note that the produced [QQ]s pair is initially almost 
pointlike, with the transverse size of rf Q ~ l/2m Q (see Fig. 3b). The produced \[QQ\ 8 g) 
state can be thus considered as a colour dipole formed by an almost pointlike colour-octet 
QQ state and a collinear gluon. The transverse size of the \ccg) state can be estimated 
[|l5f from the characteristic virtualities of the diagram of Fig. 3b as r± ~ (2m c AQco) ^ 2 — 
0.20 — 0.25 fm. An interaction inside nuclear matter will most likely prevent this state 
from binding, at later stage, to the quarkonium - the colour octet cc, with its collinear 
gluon stripped off, will preferably produce open charm mesonsQ. 

3 Another example is provided by the scattering of quarkonium states at very high energies p0[ . 
4 Note that the J/-0 production cross section represents only a tiny part, of the order of 1%, of the total 
charm production - this means that the probability to pick up a collinear gluon for the colour-octet QQ 
state is in general very small. 
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Let us try to estimate the break-up cross section of such | [QQ]sg) state in its interaction 
with nucleons |TJ|]. The transverse size of the |[cc]8<?) state estimated above is roughly the 
same as the size of J/ip. The |[cc] 8 g) state however is not bound, so, contrary to the case 
of J/ip, we do not expect any threshold suppression of the break-up cross section. We 
can therefore estimate the | [ccjgf?) break-up cross section rescaling the value of the J/ip 
break-up cross section at high energy (where the threshold suppression does not affect the 
cross section) by the colour factor 9/4, arising from the difference between the couplings 
of colour dipoles formed by the triplet and octet charges. At the energy range relevant for 
the fixed target experiments, the J/ip break-up cross section evaluated in the formalism 
of section 2.1 is <Jj/^n — 2.5 — 3 mb. We therefore get a(c Cg )N — 6 — 7 mb as an estimate 
of the | [cc] 8<7) absorption cross section. The analogous estimate for bottomonium states 
yields cr^^ ~ 1.5 — 2 mb. These estimates are admittedly rough; they show, however, 
that the nuclear attenuation of quarkonium production is in general quite strong and, in 
the first approximation, is universal for various quarkonium states. 



3. QUARKONIUM AS A PROBE OF DECONFINED MATTER 

We have shown in the previous section that the absorption cross sections of tightly 
bound quarkonium states at low energies are very small due to the softness of gluon 
fields confined inside light hadrons; this protects J/ip in a thermal hadron gas at all 
meaningful temperatures (T < 0.3 GeV) |PH|- ||14|| . On the other hand, the distribution 
of gluons in a deconfined medium is directly thermal, so that the deconfined gluons are 
hard, with the average momentum of (p 9 )deconf = 3T. An immediate consequence of 



deconfinement is thus a considerable hardening of the gluon momentum distribution [10 



IT]. Hard deconfined gluons can easily break up the J/i/j; the cross section of this "gluo- 



effect" is given by 

2tt /32\ 2 /mA 1 / 2 1 (k/e -l) 3/2 

^ = t (y) U (4) 

where k is the momentum of the gluon incident on a stationary quarkonium with binding 
energy e . We thus see qualitatively how a deconfinement test can be carried out. If we 
put a J/ip into matter at a temperature T = 0.2 GeV, then the J/ip will survive if the 
matter is confined, and will disappear if the matter is deconfined, since in the latter case 
the gluons will be hard enough to break it up. 

The latter part of this statement is in accordance with the original prediction that the 
formation of a QGP should lead to a J/ip suppression []T|j3T[] . There it was argued that in 
a QGP, colour screening would prevent any resonance binding between the perturbatively 
produced c and c, allowing the heavy quarks to separate. At the hadronization point 
of the medium, they would then be too far apart to bind to a J/ip and would therefore 
form a D and a D. Our picture complements this argument by the conclusion that 
additional suppression of physical J/ ip in dense matter will occur if and only if there is 
deconfinement. 

The dissociation of J/ip (or T) in both pictures is a consequence of the interaction with 
strong gluon fields present in deconfined matter. There is, however, a difference between 
the two mechanisms: the static screening picture takes into account the effect of deconfined 
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fields on the binding potential acting between the heavy quarks, but neglects the energy- 
momentum transfer between the J /if) and the heat bath. The dynamical "gluo-effect" 
picture of quarkonium suppression, on the other hand, emphasizes the role of the energy- 
momentum transfer from deconfined gluons to the J/ if) , but neglects the screening of the 
binding potential. Both pictures are expected to describe the physics of J/ if) suppression in 
their respective domains of applicability; they should emerge as two limits in one unified 
microscopic approach, that still has to be developed. The parameter that is relevant 
in this problem is X(T) = AE(T)/T, where the binding energy of quarkonium AE 
depends on the temperature of the system T because of the Debye screening. In the weak 
coupling limit of X <C 1, the binding energy is negligible compared to the temperature, 
and the quarkonium will simply fall apart with the rate R = 4/ ' L{T '/-kMq) 1 ^ 2 (L is 
quarkonium size), which is the classical high temperature limit of thermal activation rate 
jT2"| . In the strong coupling limit of J > 1, on the other hand, the system is tightly 



bound, and the binding energy threshold has to be overcome by the absorption of hard 
gluons from the heat bath. The rate of dissociation in this case should be computed 
from the thermal average of the gluon- quarkonium cross section (|4|). The actual value of 
X at different temperatures depends, of course, on the detailed dynamics of screening; 
lattice calculations can be of significant help here, fixing the temperature dependence of 
quarkonium mass. 

It is important to note that the dynamical "gluo-effect" approach to J /if} suppression 
does not require a thermal equilibrium of the gluon fields, so that it will remain applicable 
even in deconfined pre-equilibrium stages. Quarkonium interactions in an equilibrating 
parton gas were considered in ref. 



4. PHENOMENOLOGY OF QUARKONIUM PRODUCTION 
IN NUCLEAR COLLISIONS 4 

4.1. p — A collisions 

According to our discussion in section 2.3, in the presently accessible kinematic region 
of J /if) production by p — A collisions {xp > 0), the target nucleus sees only the passage of 
the pre-resonance state; physical charmonium states are formed outside the nucleus. The 
size of the pre-resonance state is determined by the charmed quark mass and confinement 
scale and is therefore the same for J /if) and if)' . The nuclear attenuation of J /if) and 
if)' production in p — A collisions should thus be universal. Indeed, the J /if) and if)' pro- 
duction in pA collisions shows to the same A-dependence. Fitting the available data on 
the if)' /(J /if)) ratio |4|] to the form A a leads to 

a = 0.0 ±0.02, 95%C.L.; (5) 

this rules out variations of more than 10 % between pp and pU collisions. The suppres- 
sion of J / if) production in p — A collisions should thus be understood as pre-resonance 
absorption in normal nuclear matter. This accounts naturally for the equal suppression 
observed for the two states, which would be impossible for physical resonances of such 
different sizes. 



5 This section is based on the work [33] 
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We shall now determine the pre-resonance absorption cross section from the NA38/51 
p — A data at incident proton beam energies of 200 and 450 GeV Q . In Glauber theory, 
the survival probability for a J/ip produced in a p — A collision is given by 



qGI 
D pA 



Act, 



piV— >V 



J d 2 b dzp A (b, z) exp - 1) J dz'p A (b, z')a abs ^ 



(6) 



Here p A is the nuclear density distribution, for which we take the standard three-parameter 
Woods-Saxon form with parameters as tabulated in Ref. fl34 |; it is normalized to unity, 
with J d 2 bdzp A (b, z) = 1. The suppression is thus fully determined by the absorption 
cross section cr a b s in nuclear matter. From the NA38/51 data we obtain the best fit for 



a abs = 6.3 ± 0.6 mb, 95% C.L.; 



(7) 



the corresponding survival probabilities are plotted in Fig. 4. The agreement is seen to be 
excellent in all cases. The value (]?]) is consistent with the theoretical estimates of section 
2.3, which suggest for the absorption cross section of the cc — g on nucleons a abs — 6 — 7 
mb \ TE \. 




A 

Figure 4. J/ip suppression in pA collisions; the NA38/50 data (black points) are compared 
to the Glauber theory calculations (grey points) with a a b s = 6.3 ± 0.6 mb. 



We have carried out the same analysis for the 800 GeV E772 data (see for a review); 
here the cross section is slightly larger: a abs = 7.4 ±0.7 mb , but within errors compatible 
with the value (^) obtained from the NA38/51 data. A slow increase of the absorption 
cross section with energy can be attributed to the growth of the gluon structure function 
towards smaller x (the same effect is responsible for the increase of the J/ip absorption 
cross section in the relevant energy range, see Fig. 2). 

We thus conclude that J/ip and ip' production in pA collisions is quantitatively well 
described by absorption of a pre-resonance charmonium state in nuclear matter, with the 
absorption cross section for both states in the energy range of SPS experiments given by 
Eq. (0). We now extend this description to nuclear collisions. 
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4.2. S-U Collisions 

In nucleus-nucleus collisions, charmonium production can be measured as function of 
the centrality of the collision, and hence we have to calculate the J/i/j survival probability 
at fixed impact parameter b. It is given by 



dSf B (b) 1 



d 2 b AB &nn- 



d<JAB- 

d 2 b 



d 2 sdzdz' pa{s, z)pb{P — s, z')Sa(z, s)Sb(z' , s), 



where Sa(z,s) = exp { — (A — 1) J 2 °° dzA Pa(s, za) <?abs}, and analogously for Sb(z',s). 
Here s specifies the position of the production point in a plane orthogonal to the collision 
axis, while z and z' give the position of this point within nucleus A and within nucleus B, 
respectively. The nuclear density distributions pa and p B are defined as above. To obtain 
normalized survival probability at fixed impact parameter b, we have to divide [dS Gl /d 2 b] 
by[dS Gl (b;a abs = 0)/d 2 b}. 

Experimentally, the centrality of the collision is determined by a calorimetric mea- 
surement of the associated transverse energy Et] we thus have to establish and test a 
correpondence between impact parameter b and transverse energy Et- This correlation 



can be expressed in terms of the number of "wounded" nucleons ||35|| . Each wounded 
nucleon contributes on the average an amount q to the overall transverse energy produced 
in the collision, so we have the relation 

E T (b) = qN w (b) (9) 

between the average number N w of nucleons wounded in a collision at fixed impact pa- 
rameter b and the associated average transverse energy Et produced in that collision. In 
the analysis of specific experimental results, the proportionality factor q depends on the 
details of the detector, in particular on the rapidity and transverse momentum range in 
which the produced secondaries are measured. 

The average number of wounded nucleons in an AB collision at impact parameter b is 
given by 

N% B (b) = f d 2 sn w AB (b lS )= AJd 2 sT A (s){l-[l-cr N T B (s-b)] B } + 

+ BJd 2 s T B (s- b) {l - [1 - cy N T A m A } • (10) 

Here <tn ~ 30 mb denotes the inelastic production cross section, and Ta(s) = J dz pa(z, s) 
the nuclear profile function; the s*-integration runs again over a plane orthogonal to the 
collision axis. The distribution flTU| ) is normalized in the following way: 

NJ B = — [ d 2 bN% B (b) = — (Aa B + Ba A ). (11) 

&AB J &AB 

Since there are fluctuations in the number of wounded nucleons and in the transverse 
energy of the secondaries that each wounded nucleon produces, there will be corresponding 
fluctuations in the relation between Et and b. We assume the dispersion D in the produced 
transeverse energy to be proportional to y/N w , D 2 = aErib), with a universal physical 
parameter a to be determined from pA or AB collisions. We choose the Et — b correlation 
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function Pab(-Et,^) as a conventional (see, e.g., [j36|| ) Gaussian distribution around the 
central value (^[) with dispersion D; it is normalized at fixed b: J dEj> PAB{Er,b) = 1. 
We have checked that both minimum bias [37] and Drell-Yan associated @,[17 



transverse energy spectra are very well reproduced in the approach outlined above. With 
the relation between the measured transverse energy Et and the impact parameter b of 
the collision thus determined, we can now calculate the Et dependence of the charmonium 
survival probability in nuclear matter. 

We begin with J/ip production. The experimentally determined quantity is the ratio 
[do J j[^ I dEx) I {dvAEi) of J/'ip to Drell-Yan production, measured in the mass interval 2.9 < 
< 5.5 GeV. From this we obtain the survival probability at fixed Et 



__ AB 
J exp K-^T) — jj^ 



a 



AB 



ao AB / aa AB 

dE T 1 dE T 



(12) 



by normalizing the measured ratio at fixed Et by the measured integrated cross sections. 
The quantity ( |12"1) can be directly computed in the Glauber theory formalism outlined 
above. Using the value of the pre-resonance absorption cross section (|7]), determined 
from the analysis of p — A data, we have found a good agreement with the .Sr-integrated 
O — Cu, O — U and S — U data and with Et distributions measured in S — U collisions, 
as we shall shortly show. 

For if)' production, the situation changes. The data for the integrated and the differential 
survival probabilities are considerably lower than what nuclear absorption predicts, and 
the additional suppression moreover increases with increasing Et- We therefore need to 
include the effect of additional ip' suppression on J j ip production. The branching ratio for 
the reaction ip' — > J/ip is 0.57; therefore the ip/ip' ratio measured inpp andpA collisions 
implies that 8±2 % of the observed J/ip's are due to ip' decay. Since the ip' is suppressed 
in S — U collisions, the corresponding fraction of the observed J/i/j's must be suppressed 
as well. This correction reduces theoretical predictions on the average by ~ 5%. We 
show in Fig. 5 (left) the resulting corrected theoretical Et dependence of J/ip survival 
probabilities. The agreement between the data and predictions is seen to be excellent. 

I wish to stress that we do not need to invoke any additional sources of direct J/ijj 
suppression (apart from the nuclear absorption of pre-resonance charmonium state) to 
describe the data. On the other hand, the additional ip' suppression found in S — U 
collisions clearly indicates the presence of produced matter at the stage when charmonium 
states are formed. The agreement of our Glauber calculations with the measured J/ip 
survival probabilities shows, however, that this matter cannot break up J/i/j states. This 
can be explained by the smallness of J/ip dissociation rate in confined hadronic gas |10 



12] advocated in sections 2 and 3. 



4.3. Pb-Pb Collisions 

We can further check our approach in Pb — Pb collisions, since the NA50 experiment 



17] is equipped with a zero degree calorimeter (ZDC), which determines at each Et the 



associated number of projectile spectators - those projectile nucleons which reach the ZDC 
with their full initial energy Ei n = 158 GeV/c. This additional information is important, 
since it uniquely identifies the peripherality of the collision. Denoting the projectile as A, 



11 



. 8 

Sab o.i 



. 6 
. 5 
0.4 
0.3 
0.2 
. 1 




40 



80 



. 8r 



0.7 



. 6 



U . 3 



0.5 
. 4- 
. 3 : 
.2 
. 1; 
0- 



40 



80 120 



Ex 



Figure 5. J/ip suppression in S — U (left) and Pb — Pb (right) collisions; the NA38/50 
data |§] (black points) are compared to the Glauber theory predictions |33[] (grey points) 
with <r abs = 6.3 ± 0.6 mb. 



the number of projectile spectators is evidently A — N£, with N£ of the A nucleons in the 
projectile wounded. We thus have E Z dc — (A ~ N^)E in ; using E T = q N w , we predict 
p3| the Et — Ezdc correlation which agrees very well with the measured one |I7|], 

We are now ready to address the J / ip production in the NA50 experiment [|l^] . The Pb— 
Pb results, plotted as a function of the average path L of J/ip in nuclear matter, clearly 
show strong additional suppression beyond the expected on the basis of ~ exp(—poa a b s L) 
dependence |l7j . The conclusion on the "anomalous" nature of this suppression, however, 
crucially depends on the magnitude of L, assigned to the Pb — Pb points. We would 
like therefore first to check the L assignment of the NA50 Collaboration in our approach, 
which directly gives the J ftp survival probability at a given Et- The results are presented 
in Fig. 5 (right). One can see that while the lowest Et point is still marginally consistent 
with Glauber theory, the suppression observed at higher Et indeed goes significantly 
beyond expected. Equating our calculated survival probability to the form used by the 
NA50, Sgi = exp(—poa a b s L), we confirm the NA50 L assignments flOf! - Since, as we have 
shown, the Glauber theory approach has been extremely successful in reproducing the 
bulk of J I ip production data in p — A and A — B collisions, the suppression observed in 
Pb — Pb indeed can be called "anomalous" . 



5. IS THE QUARK-GLUON PLASMA DISCOVERED? 



Before we address this provocative question, posed at this Conference also by J. -P. Blaizot 
p9| and C.-Y. Wong [|l|, let us consider the possible differences between the collision 
dynamics in S — U and Pb — Pb systems. The success of Glauber theory in describing 
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the J I if) suppression in S — U collisions and its failure in Pb — Pb points to a difference 
in the properties of matter produced in these two reactions. We shall try to describe 
this difference in terms of two variables, one of which characterizes the energy density of 
produced matter, and the other its degree of equilibration. 

In Glauber theory, the initial energy density achieved in the collision is proportional to 
the density of wounded nucleons n w (see Eq. fllDp) in the transverse plane. The average 
energy densities achieved in S — U and Pb — Pb collisions are almost identical; at first 
glance this suggests that the matter seen by produced J/ip's should be the same in both 
cases. This is not so, however, for two reasons. First, the profile of the energy density 
in the two systems is different: central Pb — Pb collisions produce a "hot core", inside 
which the energy density is higher than the highest one attainable in S — U system by 
about 25%. Second, the J/ift's are produced mostly in this central region (see Eq. (§)), 
and thus feel the matter which is hotter than average. These two effects combined lead 
to significant difference in the energy densities of matter seen by J/^'s in S — U and 
Pb — Pb collisions. This is illustrated in Fig. 6, where we plot the ratio of experimental 
J /if} suppression to the Glauber theory predictions versus the average density of matter 
seen by J/if> (to calculate this latter quantity, we convolute the density distribution with 
the J /if) production profile). Apart from being more dense, the matter seen by J /if} is also 
likely to be more thermalized. We can quantify this statement introducing the variable 

« = ^i, (13) 

w 

where v is the number of inelastic NN collisions and w is the number of wounded nucleons 
these collisions produce. The value of k tells how many times, on the average, each 
wounded nucleon was hit. In a pp collision, v — 1 and w = 2, so that k — 1. In 
pA collisions, the number of wounded nucleons in the target is equal to the number of 
collisions P5[| , which, after taking into account the wounded projectile nucleon, again 
yields k = 1. In nucleus-nucleus collisions, however, the value of k can exceed unity 
because the nucleons once wounded can collide again and again. These collisions can 
break down the independence of fragmentation of wounded nucleons and provide initial 
conditions for the onset of collective behaviour in the system. Indeed, when k > 1, 
the partons from different wounded nucleons interact, which is a necessary initial stage 
for producing deconfined matter - at large k, partons can no longer be attributed to a 
particular pair of wounded nucleons and overlap in the transverse plane. It is evident 
that in nuclear collisions k grows with atomic number and/or energy (since the number 
of collisions depends on the inelastic cross section). In central S — U collisions at SPS 
energy, the Glauber theory calculation yields R$u — 1-7, whereas for a central Pb — Pb 
collision we find RpbPb — 2.4. The Pb — Pb value thus is as far from the S — U one as 
the S — U is from pA. Central Pb — Pb collisions therefore are likely to produce not only 
more dense, but also more thermalized matter. 

It is important to note that the value of k and its variation with centrality can be 
determined in a model-independent way directly from the experimental data. Indeed, 
the number of collisions v is proportional to the number of produced Drell-Yan pairs, and 
the number of wounded nucleons to the produced transverse energy (see Eq. (||)), so at a 
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Figure 6. The ratio of experimental J/ip suppression in S — U [ffl (crosses) and Pb — Pb 
|L7| (circles) collisions to the Glauber theory predictions versus the average density of 
wounded nucleons seen by J/ip (the latter quantity is proportional to the energy density). 
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(14) 



where Ndy{Et) is the number of Drell-Yan pairs associated with a given Et bin of the 
width AEt- 

We thus conclude that the matter seen by J /if) in central Pb — Pb collisions is more 
dense and more thermalized, so the occurence of new phenomena at least cannot be 
excluded a priori If the observed J/ip suppression were indeed to be interpreted as a 
signal of deconfmement phase transition, apart from being "anomalous", it also has to 
exhibit a threshold behaviour. This feature seems to be present in the data (see Fig. 6): a 
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slight increase in the energy density induces dramatic deviation from the trend established 
previously. Moreover, if one considers the profile of the density n w and k in the transverse 
plane and assumes that all J/^'s produced in the region where n w > n™y x , k > k™{j x 
are dissociated, the resulting suppression falls below the data points by only 10 — 15% 
(similar analysis has been presented at this Conference by J. -P. Blaizot |[39|). This shows 
that the observed suppression is almost as strong as we could possibly accomodate! 

Can we still find a "conventional" explanation of the NA50 effect? It is of course too 
early to try to answer this question; let us therefore limit ourselves to some preliminary 
observations. Since Glauber calculations show that nuclear absorption of pre-resonance 
state cannot explain the Pb — Pb data, a conventional explanation has to invoke addi- 
tional suppression of J/^'s in the produced confined matter. This suppression would be 
characterized by a smooth (basically, exponential) dependence on the density of produced 
matter. Looking at Fig. 6, one realizes that it would not be easy to fit the data by this 
kind of a smooth dependence - the fits would most likely overestimate the slope of the 
S — U data and underestimate it for the Pb — Pb points. The physical reason for this is 
transparent - if hadronic comovers do not induce an additional J / ip suppression in S — U 
collisions, it is difficult to make them effective in Pb — Pb system. One may try to as- 
sume a larger density variation with E? (sometimes the density is assumed to be directly 
proportional to Et IfHI )■ Indeed, a large variation of density is possible in very central 
collisions due to fluctuations in the number of produced hadrons. However the measured 
Et — Ezdc correlation discussed in section 4.3 shows that in the presented Pb — Pb data 
the variation of Et results from the variation of the collision centrality. Even the highest 
Et point of NA50 is not entirely in the fluctuation domain, and corresponds to the mean 
impact parameter of b ~ 2 fm. In this regime, the Et measured by the NA50 varies by 
more than four times, but this leads only to ~ 30% variation of the initial energy density 
(see Figs. 5 and 6). Additional constraint on a conventional scenario is imposed by the 
ip'/tp ratio: if the density of comovers in Pb — Pb were much higher than the density in 
S — U, this would imply a smaller ip/i/j' ratio in the former case - the prediction that 
would bring us in conflict with the data [IT]]. A possible way to describe the Pb — Pb 
points conventionally would be to decrease artificially the value of the nuclear absorption 
cross section a a bs, leaving thus room for comover absorption already in S — U collisions 
(the need for comover effects to explain the S — U data was advocated at this Conference 
by S. Gavin [[11]). This, however, would contradict to the pA data, that fix the value of 
a a bs rather precisely (see section 4.1). Of course, it remains to be seen if a convincing 
conventional explanation can eventually be found. 

To summarize, the NA50 Pb — Pb results indeed seem to suggest that a new mechanism 
of J/ip suppression sets in at higher energy densities. The observed effect can be considered 
as a strong evidence of some kind of deconfinement in nuclear collisions. What can we do 
to turn this evidence into a proof, or to discard it? 



6. WHAT HAS TO BE DONE NEXT? 

New precision data on quarkonium production coming from CERN SPS, Fermilab, 
HERA and elsewhere allow us today to get rid of many uncertainties inherent to the 
analyses of J/tp suppression over the years. A coherent picture, providing a good descrip- 
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tion of the bulk of existing pp, pA and AB data, has started to emerge - this makes us 
ready to recognize and study unusual phenomena. It is therefore particularly important 
to learn more about the onset of anomalous behaviour of J/ip suppression seen by the 
NA50 Collaboration. More statistics and more data points, both in the transition regime 
of small Et and in the fluctuation region of the highest E?, are needed to establish the 
threshold behaviour suggested by the present data. An important information would be 



also provided by the J/tjj transverse momentum distributions |j42 |. 

Direct measurement of the low-energy J/ip absorption cross section in the proposed 
inverse kinematics experiment |10|,[1^] has become possible with the advent of Ph beam 



at CERN SPS. This experiment would allow us to directly constrain the J ftp absorption 
possible in a hadronic medium, providing important additional check of the deconfmement 
transition as the cause of "anomalous" J/ip suppression. 

Heavy quarkonium represents a rare example of a strongly interacting system that is 
simple enough to be systematically analyzed by the current theoretical methods. It has 
already proved to be extremely useful for understanding the properties of QCD and its 
ground state - the vacuum. I believe that quarkonium will tell us much also about the 
critical behaviour of QCD matter produced by relativistic heavy ion collisions. 
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